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ABSTRACT 



We present an analysis of infrared (IR) echelle spectra of five stars in the TW 
Hydrae Association (TWA). We model the Zeeman broadening in four magnetic- 
sensitive Ti I lines near 2.2 /im and measure the value of the photospheric mag- 
netic field averaged over the surface of each star. To ensure that other broadening 
mechanisms are properly taken into account, we also inspect several magnetically 
insensitive CO lines near 2.3 /im and find no excess broadening above that pro- 
duced by stellar rotation and instrumental broadening, providing confidence in 
the magnetic interpretation of the width of the Ti I lines. We then utilize our 
results to test the relationship between stellar magnetic flux and X-ray properties 
and compare the measured fields with equipartition field values. Finally, we use 
our results and recent results on a large sample of stars in Taurus to discuss the 
potential evolution of magnetic field properties between the age of Taurus (~2 
Myrs) and the age of TWA (~10 Myrs). We find that the average stellar field 
strength increases with age; however, the total unsigned magnetic flux decreases 
as the stars contract onto the main-sequence. 

Subject headings: infrared: stars — stars: magnetic fields — stars: pre-main 
sequence — stars: individual (TWA 5A, TWA 7, TWA 8A, TWA 9A, TWA 9B) 
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Introduction 



T Tauri stars (TTSs) are young, low-mass stars, which are newly formed out of their 
natal molecular clouds and are often associated with OB associations and/or nebulosity. 
The location of TTSs above the main sequence in the H-R diagram and their lithium 
abundance also reveal their youth. TTSs display substa ntial photometric and s p ectro scopic 



variabili t y, and their properties have b een r eviewed by 



Bertout 



fll989h 



Basri fc Bertout (il993l ) and 



Appenze 



Menard fc Bertout 



ler fc Mundt 



fll989h 



19991 ). Depending on 



whether they are actively accreting, TTSs are divided into two categories. Classical T Tauri 
stars (CTTSs) are still surrounded by dusty circumstellar disks with detectable non-stellar 
emission in the near-infrared (near-IR). Accretion from the disks onto the stars produces 
many observable features such as strong yet variable Ha emission, often displaying P-Cygni 
or inverse P-Cygni profile shapes. Naked T Tauri stars (NTTSs), on the other hand, show 
little or no sign of accretion and do not appear to have close dusty disks around them. It is 
generally believed that stars evolve through a CTTS phase to become a NTTS on the way 
to the main sequence. 

One essential piece of the puzzle to our understanding of TTSs is the stellar magnetic 
field. For NTTSs, it is believed that dynamo generated magnetic fields are responsible 
for most of their peculiarities, e.g., variable X-ray emission, flares, and cool spots. The 
evolution of CTTSs, especially the interaction between the stars and their circumstellar 
disks, is believed t o be directly controlled by stellar magnetic fields (see review by 



Bouvier et al 



1991 



20071). Various magnetosp 



Collier Cameron fc Campbell 



leric accretion models ( 



1993 



Shu et al. 



1994 



Camenzind 



199G 



Paatz fc Camenzind 



Konig] 



19961) 



generally agree that the disks of CTTSs are truncated by stellar magnetic fields near the 
corotation radius and the accreting material fiows along the field lines toward the polar 
regions of the central star. The disks of CTTSs are the birthplace for planets and the 
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radius at which the magnetic field truncates the inner disk may determine where giant 



planets are halted on their inward migration in the disk (ILin et al 



1996 



Eisner et al. 



20051). 



Magnetic fields probably a lso play a key role in driving the stellar winds/jets seen in young 



stars (e.g., 



Shu et al. 



19991 ). Therefore, a full understanding of the early evolution of young 
stars, which will lead to better understanding of both star and planet formation, requires 
detailed theoretical consideration and empirical measurements of the stellar magnetic field 
properties. 

Direct detection of stellar ma gnetic fields is diff icult, with measurements starting 



to become plentiful only recently (jJohns-KruU 



20071 ). One 



first direct detection of magnetic fields on a TTS, Tap 35 (IBasri et al. 



technique, which led to the 



other stars ( IGuenther et al 



19921 ) and later on 



19991 ). is based on an equivalent width increase in moderately 
saturated magnetically sensitive lines due to Zeeman splitting. Unfortunately, this method 
is particularly sensitive to the choice of effective temperature and gravity used when 
analyzing the stellar line measurements. Another method of detecting stellar magnetic ffelds 
is through spectropolarimetry, which looks for net circular polarization in Zeeman-sensitive 
spectral lines. Observed along the direction that is parallel to the magnetic ffeld, the 
Zeeman a-components are circularly polarized, with the components of opposite helicity 
split to either side of the nominal line wavelength. If the magnetic topology on the stars 
is somewhat organized, there should be a measurable component along the line of sight. 
Reliable detections of circular polarization in TTS photospheric abs oprtion lines are ra re, 



with most recent studies 


1999a 




1 

Daou et al.1 


20061; 1 



weak, detections have been made (lYang et al. 



20041) . More recently, signiffca nt, but similarly 



2007 



Donati et al 



strong detection on BP Tau has recently been published ( IDonati et al. 



2007) 



ana 



one surprismg 



20081 ). On the other 



hand, net polarization is detected in the He I 5876 A emission line, ffrst discovered by 



Johns-KruU et al 



1999al ) on BP Tau. This He I line is believed to form in the postshock 
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region (iHartmann et al.lll994j : Edwards et al.lll994j ) where material accretes onto the star 
and the detection on BP Tau yielded a net longitudinal field of 2.46 ± 0.12 kG in the 
line formation region. Circular polariz ation in the narrow component of the He I line has 



now been o 



Yang et al 



jserved in several CTTSs ( Va 



2007 



Donati et al. 



2007 



enti fc Johns-KruU 



2004 



Symington et al. 



2005 



20081 ). These observations suggest that accretion onto 



CTTSs is indeed controlled by a strong stellar magnetic field. 

While spectropolarimetry can provide clear detections of stellar magnetic fields, it is 
only sensitive to the net field resulting from the predominance of one field polarity over the 
other. As a result, it can miss the majority of the magnetic fiux present on a star. The 



the Zeeman 


broadening of spectral lines in unpolarized light ( 


Robinson 


1980; 


S 


Valenti et al. 


1995; 


Johns-KruU & Valenti 


1996; 


Johns-KruU et al. 


19991 


) ) . For 



Saar 



1988 



Zeeman component, the splitting resulting from the magnetic field is 



AA 



:y gB = A.67 X 10~'y gB mA 



where g is the Lande ^f-factor of the transition, B is the strength of the ma gnetic field in kG 



Johns-KruU et al. 



(jl999bl . hereafter Paper 



and A is the wavelength of the transition in A. 
I) detected Zeeman broadening of the Ti I line at 2.2233 fim on the CTTS BP Tau and 
obtained a field strength of i? = 2.6 ± 0.3 k G, averaged over th e entir e surface of the star. 



From analysis of four Ti I lines near 2.2 fim, 



Johns-KruU et al 



(120041 . hereafter Paper II) 



also detected a strong magnetic field of 2.5 ± 0.2 kG on the NTTS Hubble 4. In addition, 
these authors examined several CO lines near 2.3 fim, which are magnetically insensitive, 
and found no excess broadening above that due to rotational and instrumental effects. T his 



technique has now been used to rn easure the fields of 16 TTSs (1 Johns-KruU et al. 



Yang et al 



2004 



2005 



Johns-KruU 



20071 ) 



In the largest study to date of fields on TTSs, 



Johns-KruUl (120071 ) used a sample 
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of 14 stars to look for correlations between the measured magnetic field properties and 
stellar propertie s that might be ini portant for dynamo action. No significant correlations 



were found and 



Johns-Krulll (120071 ) speculated that the fields seen in these young stars 



may be entrained interstellar fields left over from the star formation process (e.g., 



1987 



Moss 



Tayler 



20031 ). To supply more observational constraints and further investigate the 
magnetic properties of TTSs, we present here Zeeman broadening measurements of five 
NTTSs in the TW Hvd rae association fTWA) : TWA 5 A, TWA 7, TWA 8 A, TWA 9A and 



TWA 9B. At ~ 50 pc flWichmann et al 



formation. 



Kastner et al 



19981 ). TWA is the nearest region of recent star 



(119971 ) first concluded that a few T Tauri stars in the apparent 
neighborhood of TW Hya are indeed physically close to each other, based on similar 
X-ray emission and lithium abundance, and proposed the existence of the TW Hydrae 
Association. Over 20 stars have been proposed as members of TWA (though only a few of 



all proposed TWA members have trigonometric par a. 



of TWA members can be found in 
at about 10 Myrs (see discussion in 



Zuckerman et al 



lax c onfirmation), and a compilation 



(2004) 



Barrado Y Navascue; 



The age for TWA is estimated 



20061 ). and it is intriguing that 



these pre-main-sequence stars are not associa ted with 



the TWA is still under investigation (e.g., see 



any ra olecular cloud. The origin of 



Makarov 



20071). 



As describe d earlier, the mag netic field properties of 



investisgated by 



Setiawan et al 



Yang et al 



(I2OO5I . hereafter. Paper III) and 



W Hva (TWA 1) have been 



Yang et al. 



(I2OO8I ) report a massive planet very close to TW Hya, and 



(120071). Interstinglv . 



Huelamo et al. 



( 120081 ) argue that the observed radial velocity variations of TW Hya are best explained 
by a cool spot instead of a hot jupiter. Here, we study the magnetic fields of 5 additional 
members of TWA, utilizing infrared (IR) spectroscopy to measure Zeeman broadening in 4 
Ti I lines in the K band. The observations and data reduction are presented in § 2. In § 3 
we describe our data analysis technique and results, and it is followed by a discussion of the 
results in 5 4. 
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2. Observations and Data Reduction 

We obtained K-band IR spectra of the TWA stars at the NASA Infrared Telescope 
Facihty (IRTF) in January 2000. Observa tions were made with the CSHELL spectrometer 



( iTokunaga et al 



1990 



Greene et al 



19931 ) using a 0'/5 slit. With this slit width CSHELL 
normally delivers a spectral resolving power of R = A/AA > 40,000. At the time our 
observations were taken, the resolution achieved with the O'.'S slit was degraded due to 
problems within the instrument. Soon after our observing run, CSHELL was serviced to 
restore its spectral resolution to its nominal value. Our observations achieve an actual 
resolution of ~ 24, 000, corresponding to a FWHM of ~ 4.5 pixels on the 256 x 256 InSb 
array detector. A continuously variable filter (CVF) isolated individual orders of the echelle 
grating, and a ~ 0.0057 /xm portion of the spectrum was recorded in each of 3 wavelength 
settings. The first two settings (2.2218 and 2.2291 /im) each contain 2 magnetically sensitive 
Ti I lines. The third setting (2.3125 /im) contains 9 strong, magnetically insensitive CO 
lines. Each star was observed at two positions along the slit separated by 10". Custom 
IDL software described in Paper I was used to reduce the K-band spectra. The reduction 
includes removal of the detector bias, dark current, and night sky emission. Flat-fielding, 
cosmic ray removal, and optimal extraction of the stellar spectrum are also included. 
Wavelength calibration for each setting is based on a 3rd order polynomial fitted to nX for 
several lamp emission lines observed by changing the CVF while keeping the grating position 
fixed. Dividing by the order number, n, for each setting yields the actual wavelength scale 
for the setting in question. In each wavelength setting, we also observed a rapidly rotating 
early star to serve as a telluric standard. All spectra shown here have been corrected for 
telluric absorption. The IR observations are summarized in Table [H 
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3. Analysis and Results 

The analysis technique we employ here follows that in Papers I-III where additional 
details can be found. We model the profiles of four Ti I absorption lines in the K band 
and measure the mean magnetic field (averaged over the entire stellar surface) for the 
five TWA stars. Diagnostics in the IR are generally more sensitive to Zeeman broadening 
than lines in the optical thanks to the dependence of magnetic splitting (see eq. [1]) 
compared with the dependence of Doppler broadening mechanisms. Using estimates of 
the key stellar parameters of our objects, namely, Teg, logs', [M/H] and vsini, we first 
construct model atmospheres and synthesize the Ti I line profiles without the presence of 
magnetic fields. The observed magnetically insensitive CO lines are also synthesized to 
show that instrumental and rotational broadening are properly taken into account and the 
non-magnetic spectrum is accurately predicted. There are no free parameters involved in 
computing the non-magnetic line profiles. The excess broadening seen in the Ti I lines 
relative to the non-magnetic model is interpreted as Zeeman broadening. Fitting the 
magnetically broadened Ti I line profiles then yields the average magnetic field strength on 
the stellar surface. 

Paper III presented extensive tests of this analysis technique to quantify possible 
systematic errors resulting from inaccurate stellar parameters. The Monte Carlo analysis 
showed that a 200 K error in T^s or 0.5 dex error in log (7 introduces only about 10% or less 
systematic error in the derived mean field measurements. This is because we are modeling 
actual Zeeman broadening instead of some other effect such as a change in equivalent width 
which is particularly sensitive to errors in Tes and logs'. Thus, instead of performing a 
detailed spectral analysis to determine the stellar parameters as in our previous studies 
(e.g.. Papers I-III), we adopt the most accurate values available in the literature for our 
targets. Values for these non-magnetic stellar parameters are given in Table [2] and justified 
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in subsequent subsections for each star. 



The CO hues near 2.3 /xm serve as an excellent verification of our chosen stellar 
parameters. For each of o ur targets, we interpolate a grid of the "next generation" 



(NextGen) atmospheres of 



AUard fc Hauschildt 



(I1995I ) to the specific values of Tes, log (7 



and [M/H] and construct a model atmosphere appropriate for the star. We then synthesize 
the CO line profiles and convolve them with a Gaussian corresponding to the resolution of 



Goorvitchl (jl99J) for the CO X S^ 



the observed spectrum. The CO line data are taken from 
lines. No stellar parameter fitting is involved with the synthesis of the CO lines. Continuum 
normalization is the only free parameter used to match the observed and synthetic spectra. 
The shape, especially the widths, of the synthetic line profiles is solely determined by the 
stellar parameters. A good match between the widths of the synthetic CO lines and that 
of the observations, as seen in the bottom panels of Figures 1 and 3 — 6, shows that the 
dominant non- magnetic broadening mechanisms are properly taken into account. 



A polarized radiative transfer code (jPiskunov 



I999I ) is utilized to compute the line 



profiles for both the CO and Ti I lines, assuming a radial magnetic field geometry in the 
photosphere. To construct our models, we treat the stellar surface as if it is divided into a 
number of regions, where one region is field-free and the other regions are each covered by 
a single magnetic field strength. The fractional size of each region is specified by a filling 
factor and the total of the filling factors must be unity. These regions are assumed to be 
spatially well mixed over the surface of the star. We compute the spectrum for each region, 
and the spectra of all regions are multiplied by their corresponding filling factors and 



summed up to form the final model spectrum for a specific set of fie 



facto rs. We use the nonlinear least-squares technique of Marquardt (IBevington fc Robinson 



d strengths and filling 



I992I ) to solve for the best-fit combination of field strengths and filling factors as described 
below. 
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For each target, we fit three different models. The first model (hereafter, Ml) allows 
a single magnetic field strength (B) to cover a fraction of the stellar surface (/), and the 
rest of the surface is field-free. The mean stellar magnetic field, B, is just Bf. The free 
parameters of this model are B and /. The second model (hereafter, M2) divides the 
surface of the star into three regions: two regions covered with magnetic fields of different 
strengths {Bi, B2) and filling factors (/i, /2), and one region with no magnetic field. The 
mean magnetic field on the stellar surface for M2 is, B — Bifi + i^2/2- The free parameters 
of this model are Bi, /i, B2 and /2. The third model (hereafter, M3) allows three magnetic 
regions on the stellar surface to have specific field strengths of 2 kG, 4 kG and 6 kG, 

respectively. The rest of the surface is field-free. Only the filling factors are solved for in 

3 

this model. The mean field is then B — B^fi. 

i=l 

As discussed in Paper I-III, the Ml model provides a much better fit to the observations 
than models with no magnetic broadening, but the model Ml profile is not as smooth and 
broad as the observed line profile, and is incapable of supplying good fits for the cores and 
the wings of the lines simultaneously. The magnetic structure on the surface of TTSs is 
unlikely to be so simple, and this is the motivation to add more magnetic components in 
model M2 and M3, which provides improved smoother fits to the observed fine profiles. 
Due to lower spectral resolution (~ 24,000) than that in Paper I-III (~ 36,000), the 
need for more magnetic components is less urgent. For example, likely due to the lower 
spectral resolution achieved here, the data do not clearly require a model with 3 magnetic 
components as in M3; however, we fit such a model to be consistent with previous work in 
order to avoid any systematic differences when comparing the field measurements of TWA 
stars with previous results. Below, we discuss the results for each of our targets in detail. 
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3.1. TWA 5A 



TWA 5 is am ong the five original TWA members reported by iKastner et al.l (119971 ) 



Webb et al 



(119991 ) discovered a M8. 5 brown dwarf co n ipanio n, TWA 5B, ~ 2" away from 



tfie mucfi brighter TWA 5A (Ml. 5). 



Macintosh et al 



(I2OOII ). using adaptive optics (AO) 



imaging, found that 
~ 0'.'06 at that time. 



WA 5 A itself is also a close binary system with a separation of 



Konopacky et al. 



(120071 ) monitored the TWA 5Aab system for five 



years. Their speckle and AO observations yield an orbital solution with a semi-major axis 
of 0'.'066 ± O'.'OOS and a period of 5.94 ± 0.09 yrs. These authors also find a K band flux 
ratio of ~ 1.25, so we are most likely seeing substantial light from both stars in our spectra. 



It is also suspected that there is possibly another s pectroscopic 



comp onent in the 



2OO3I). We do notice 



TWA 5Aab pair, causing large radial velocity variations (iTorres et al.l 
peculiarities in our spectra of TWA 5Aab system, which is shown in Figure [H The three 
wavelength settings shown here were observed on 3 consecutive nights. The two Ti I line 
profiles in the 2.2300 /xm setting are very asymmetric, and some of the CO lines show large 
wavelength shifts or asymmetries as well. With a period of about 6 years, the line profiles 
of the binary system obtained on consecutive nights should not vary so dramatically. This 
drastic change may be caused by a third component in the TWA 5Aab system in a close 
orbit with one of the two known members. 

We do not notice obvious asymmetry in the two Ti I line profiles in the 2.2228 fim 
setting, so these two lines are used for magnetic field measurement. Following 



Weintraub et al. 



(l2000l ). we adopt an effective temperature of 3700 K and derive a log (7 of 



3.7 from t heir absolute magn itude and stellar mass estimates, which assume a distance o f 



55 ± 9 pc (IWebb et al. 



I999I ). A vsini of 38 km s ^ is adopted from 



Torres et al. 



(12003h . 



The measured mean magnetic field values of TWA 5A along with that of the other four 
TWA stars, using the three different models, are listed in Table [21 The reduced values 
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for TWA 5A given in Table [2] are calculated with all four Ti I lines included so that it is 
consistent with other stars. Note that only the two lines in the 2.2228 fim region are fitted 
and the reduced values calculated with these two lines only are 1.06, 1.10 and 1.28 for 
models M1-M3, repectively. The observed and the best fit model spectra of TWA 5A are 
shown in Figure [TJ Given the uncertainties associated with the number of stars in this 
system, we report our field estimate here (which is likely some sort of average of all the 
stars present) but do not use measurements of this system in any subsequent analysis. 



3.2. TWA 7 



TWA 7 was discovered as a TWA member by I Webb et al.l (119991 ) . iNeuhauser et al. 



(I2OOOI ) detected a possible extra-solar planet candidate about 2'.'5 away from TWA 7 



using ground-based direct imaging, but they concluded that this object is more likely an 



unrelated back ground source 



observations of 



Lowrance et al. 



lis in terpretation was later confirmed by the coronagraphic 



(120051 ). An inchnation of 28° is estimated for TWA 7. 



Webb et al. 



(I1999I ) determined a spectral type of Ml 



low resolution optical spectra with spectral standards from 



or TWA 7 by comparin g their 



Montes fc Martin 



(119981 ). In our 



high resolution IR data, we notice that some molecular lines of TWA 7 in the K band are 
more prominent than that indicated by a Ml spectral type. Shown in Figure [21 we compare 
the spectrum of TWA 7 in the 2.2228 /im setting with that of several stars of different 
spectral types. From this TWA 7 appears to be of later spectral type than Ml, probably 
closer to M3. Though we do not synthsize molecular lines in this wavelength setting, the 
Ti I lines are not seriously blended with these features. The molecular lines, possibly CN, 
do not substantially affect our magnetic field measurement. In our analysis we adopt a 
TefT of 3300 K for TWA 7, appropriate for an M3 star, and a lo g q = 3.85 is derived by 



adopting the bolometric luminosity and stellar mass estimates of 



Neuhauser et al. 



(l2000h 
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he c onversion of spectral type to effective temperature in this paper is based on 



Johnson 



( 119661 ) ■ Interestingly, Joh nsons spectral 



derived for young stars by 



Luhman et al. 



type t o Teff conversions closely parallel those 



(I2OO3I ). from G through mid-M spectral types. We 



use a V smi = 9.0 km s ^, which is s lightly larger than previous measured value of 6.1 ±0.5 



km s ^ from 



Neuhauser et al. 



(120001 ). Our choices of f sini are based initially on literature 
values and in some cases they are adjusted by a small amount to better match the width 
of the observed nonmagnetic CO lines. As mentioned above, our analysis technique suffers 
~ 10% or less systematic error for a 200 K error in Tcfj or 0.5 dex error in \ogg, so the 
measured field values of TWA 7 are not subject to much change, whether its spectral type 



is Ml or M3. The best fit for TWA 7 is shown in Figure [31 



3.3. TWA 8A 



TWA 8 A was also identified by 



Webb et al 



(I1999I ) as an association member, and 



is classified as a M2 star. Its fainter companion TWA 8B is about 13" to the south 
of the primary. The membership of the companion in TWA is based primarily on 
its similar lithium abundance and radial velocity compared to other members of the 
association. T rigonometric parallax a nd pr oper motion measurements are needed for further 



confirmation. 



White k Hillenbrand! (120041 ) derive a spectral type of M3 (±0.5) for TWA 



8A, by matching line and molecular band ratios to those of dwarf spectral standards which 
are rotationally broadened. These authors also determined an upper limit of 5.5 km s~^ for 
the vsmi of TWA 8A. For our analysis, we use a T^s of 3400 K, \ogg of 4.0, and vsmi of 
4.0 km s~^. The inclination of TWA 8A is estimated to be 16°, so it is seen close to pole-on. 
The best fit for the Ti I lines as well as the predicted CO spectrum of TWA 8A are shown 
compared with the observed spectra in Figure HI 
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3.4. TWA 9A & 9B 



TWA 9 (CD -36°7429) was found by Ijensen et al.l fllQQSh to be in the physical 
vicinity of T W Hya. It is a binary system with a separation of 6". The Hipparcos satelhte 



flESA 



19971 ) measured a distance of 50.3 ± 6.0 pc for TWA 9 syste m. In the hterature, 



TWA 9 A is classified either as K5 f Webb et al. 



2004r). Previously reported f sin i values are 7 ± 3 km s ^ ( lReidll2003l ) and 11 ± 1 km s ^ 



19991) or K7 (±1) (IWhite fc Hillenbrand 



( jWhite fc Hillenbrand 



200J) for TWA 9A. We adopt T^s = 4000 K and log^ = 4.5, and use 



vsini = 9 km s ^. The inclination angle of TWA 9A is inferred to be 90°. 



For TWA 9B, reported 



( IWhite Sz Hillenbrand 



spectral types include Ml (iWebb et aL 



20041). and vsini v alues of 4 ± 3 km s ^ ( iReid 



1999h an d M3.5 (± 0.5) 



9 ± 1 km s"^ ( White fc Hillenbrand 



20031 ) as well as 



20041 ) have been published. We use Tes = 3400 K, 



logg = 4.0, and vsini = 10 km s~^. The inclination angle is estimated to be 60°. Due to 
the very poor signal-to-noise ratio we obtain on this star, the spectrum in the 2.2300 fim 
setting is excluded from our analysis. The best fit spectra along with observed spectra of 
TWA 9A and TWA 9B are shown in Figure [5] and O respectively. 



4. Discussion 

By fitting the Zeeman broadening in the magnetically sensitive Ti I lines, we measure 
the average surface magnetic field of five TWA stars. Three different models were used to 
fit the line profiles and measure the magnetic field strength with generally good agreement 
in the results between the three. Model Ml assumes only 1 magnetic component to the 
stellar atmosphere. This model provides a good match to TWA 5A in the 2.2228/im region 
(but not the 2.2300/im region). This is likely due to the fact that TWA 5A apparently has a 
relatively large vsini so that rotational broadening, in addition to the magnetic broadening. 
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strongly affects tlie line profile shapes. Because our magnetic models for TWA 5A do not 
reproduce line profiles in the 2.2300/im region, our derived magnetic properties for TWA 
5A may have large systematic errors. Generally, mean field strength increases as more 
components are added to the model, which suggests that the observed profiles have wings 
that are too broad to be fit by a single component. For the other stars, the lowest reduced 

typically occurs for model M2, though M3 is favored in one case. However, there is 
generally quite good agreement bewteen the mean field values recovered using these two 
models with the difference bewteen them averaging 15%. For consister icy with previous 



studies, we use the results from M3 when comparing to results found by 



■Tohns-Krulll mm . 



Theoretically, equipartition arguments for the Sun and cool stars suggest that the 
magnetic field strength should scale with the gas pressure in the surrounding non-magnetic 



photosphere which sets a limit for t he maximum fie. 



Spruit fc Zweibel 



1979 



Safier 



19991). 



Rajaguru et al. 



d stre ngth allowed on the star (e.g. 
(I2OO2I ) found that in the regime where 



convectively stable flux tubes exist, equipartition pressure allows a maximum magnetic field 
strength of ~ 1300 G for values of T^e and \ogg appropriate for most of our TWA stars. We 
calculate equipartition field strengths for the TWA stars by balancing magnetic pressure, 
B'^g/iSn), with the gas pressure in surrounding unmagnetized photosphere, Pg. We use the 
gas pressure at the level in the atmosphere where the local temperature is equal to the 
effective temperature in the model atmospheres used for our stars. This should be an upper 
limit for gas pressure because this level is approximately where the continuum forms, while 
the Ti I lines form over a wide range of depths above this level in the atmosphere. The 

1 /2 

equipartition field strengths, B^q = (SnPg) , of the TWA stars are listed in the last colum n 



Rajaguru et al. 



of Tab le [2l These values are in good agreement with the estimate from 
( 120021 ). Compared with the measured mean field values, they are all smaller by a factor of 
at least 1.38 and more typically ~ 2.0 (for model M2), suggesting that magnetic pressure 



dominates over gas pressure in the photospheres of the TWA stars observed here. This 
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in turn implies complete coverage of magnetic fields over surface of these stars. Magnetic 
regions with excess pressure would expand, while field free regions contract into the lower 
pressure nonmagnetic regions until equipartition is established at a lower field strength or 
until the entire star is covered with magnetic fields. As discussed in Papers I and II, our 
models do not necessarily require field-free regions. The results are essentially the same 
with or without a kG field component. As an example, we fit the observed spectra of 
TWA 8A with a model which allows on the stellar surface no field-free region and four 
magnetic regions with specific field strengths of 1 kG, 3 kG, 5 kG and 7 kG. The measured 
magnetic field value using this model is 3.4 kG with a reduced of 1.22, while model M3 
yields 3.3 kG with a reduced of 1-25 (see Table [2]). With a field-free region included in 
our current analysis, all the TWA stars have at least 76%, in most cases 80% — 86%, of 
their surfaces covered with magnetic fields. The surfaces of TTSs are quite possibly fully 
covered by strong magnetic fields. 



Pevtsov et al 



(120031 ) examine the X-ray luminosities and the total unsigned magnetic 
flux of active stars as well as regions of different activity level on the Sun. They find a 
nearly linear correlation over many decades of variation in the two quantities. We calculate 
the magnetic flux of the four TWA stars (excluding TWA 5A since it is an unresolved 
multiple system), using our magnetic field measur e ment s and the stellar radii estimated 



from bolometric luminosities given in 



Webb et al 



(119991 ) and the effective 



temperatures 



of our stars. We utilize the X-ray luminosity-magnetic fiux relationship of 



Pevtsov et al. 



( I2OO3I ) and calculate the predicted X-ray luminos ity. For conipariso n, we use the X-ray to 



Webb et al. 



(119991 ) to get the value of 



bolometric luminosity ratio, Lx/Lboi, reported by 
for each star. The predicted and measured values of for the three TWA stars (TWA 
7, TWA 8A and TWA 9A) presented in this paper and that of TW Hydrae taken from 
Paper III are plotted in Figure [71 A measured value for TWA 9B is not available in the 
literature. As shown in Figure. U\ the observed L^. values of the four TWA stars are smaller 
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th an the predicted va lues by a factor of 2.7 ± 0.7 on average. This is similar to the findings 



of 



Johns-Krulll (120071) in a study of 14 CTTSs in Taurus and Auriga. The stars from 
Johns-Krulll ( 120071 ) are also plotted in Figure U\ and the difference between the observed 
and the predicted values is a bit more obvious. The average ratio of the predicted to 
observed values is 17.7 ± 3.4 for the 14 Taurus and Auriga stars. At an age of about 10 
Myrs, the TWA stars are substantially older than the stars in the Taurus-Auriga region ( 
~2 Myrs). As a result, the TWA stars are substantially closer to the main sequence and 
this may be the reason why their X-ray luminosit ies and magnetic flux values are closer 



to the relationship found by 



Pevtsov et al. 



( 120031 ) which was based primarily on main 



sequence stars. The X-ray emission on active stars such as TTSs is thought to mainly result 



( Giidel et al. 


2003; 


Arzner et al. 


2007) 



20071 ). Magn etic fields are stressed due to convective gas 



20041 ). producing flares of various 



motions in the lower photosphere and below (iFoukall 
scales. For PMS stars that are fully covered by magnetic fields, when the magnetic field 
strength is larger than the equipartition value, it may be more difficult for the gas, at least 
in the photosphere, to push around the fields and build up magnetic stress. This would 
be the case for both the Taurus stars and the TWA stars. The ratio of the observed to 
equipartion magnetic field, Bobs/Beq, for the Taurus stars on average is greater than that 
for the TWA stars, so the X-ray production might be more prohibited for the Taurus stars 
resulting in the observed being further away from the line of equality in Figure [71 It will 
be interesting to see where TTSs in other age groups fit in this picture relating magnetic 
flux and X-ray emission. 

The origin of magnetic fields on TTSs has sparked various theoretical speculations, 
but it is still far from clear how the strong fields on these stars are produced. Due to their 
youth, most TTSs are generally thought to be fully conv ective, though mode l calculations 



have shown that TTSs might have small radiative cores (iGranzer et al 



2000). As a result. 
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a solar-like interface dynamo (e.g., 



Durney fc Latour 



1978 



Parker 



19931 ) is unlikely to be 



responsible for the fields on TTSs. One potent ial alternative mechanism is a distribut ed 
dynamo amplified by co nvective motions (e.g., 



Chabrier fc Kiiker 



Durney et al 



19931 : 



Dobler et al 



2006 



20061 ). Such turbulent dynamos are not well understood at this point. 



but exisiting models seem to only produce field strengths equal to or lower than the pressure 



equipartition values, and many such models a. 
small filling factor at the stellar surface (e.g.. 



so produce strong fields which occ upy only a 



Cattaneo 



1999 



Bercik et al. 



20051 ). Both of 



these predicitons appear not to be borne out by the observations. In addition to a possible 



dynamo origin, the magnetic fields of T TSs cou 



survives the star formation p rocess (e.g., 



from this paper and those in 



Yang et al. 



'avler 



d also resu 



1987 : 



(j2005|) and 



Moss 



; from the primordial flux that 



2003). Wit 



Johns-KruU 



1 the measurements 



(120071 ). we explore this 



issue by comparing the average magnetic field strength and magnetic flux of the stars in 
TWA and in Taurus-Auriga, two regions of different ages. The average magnetic field of 
the Taurus-Auriga stars is 2.20 ± 0.16 kG, which is smaller than that of the five TWA stars 
(TW Hya, TWA 7, TWA 8A, TWA 9A and TWA 9B), where B = 2.98 ± 0.22 kG. While 
the average field strength increases slightly from the age of Taurus to that of TWA, we see a 
sharp decrease of the average mag netic flux, inR.^B, from (7.00 ± 1.37) x 10^^ Wb for the 
Taurus-Auriga stars compared to (2.61 ± 0.69) x 10^*^ Wb for the five stars in TWA. This 



change in the magnetic flux is general 



fields on TTSs as suggested by 



y con sistent with the decay of primordial magnetic 



Mosd (120031 ) . With only a limited number of samples on 



hand, the current evidence is far from sufficient. Further observations and measurements of 
young stars in other young clusters of different ages, e.g., the Orion Nebula Cluster and the 
P Pictoris Moving Group, could shed more light on the issue of field generation in young 
stars. 
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Table 1: Journal of Observations 



Object WavoloriRth Sottinp; UT dato UT time Total Exposure Tiine(s) 

TWA 5A 2.2228 /xm 2000 Jan 08 14:08 2400 

2.2300 /xm 2000 Jan 09 14:51 2400 

2.3129 Atm 2000 Jan 07 15:03 2400 

TWA 7 2.2228 /xm 2000 Jan 08 13:22 2400 

2.2300 nm 2000 Jan 09 13:04 2400 

2.3129 Atm 2000 Jan 07 14:00 3000 

TWA 8A 2.2228 /xm 2000 Jan 08 15:01 2400 

2.2300 Atm 2000 Jan 09 13:50 2400 

2.3129 /xm 2000 Jan 11 11:43 3600 

TWA 9A 2.2228 /xm 2000 Jan 11 11:43 3600 

2.2300 Atm 2000 Jan 12 13:26 3600 

2.3129 iim 2000 Jan 11 13:24 4800 

TWA 9B 2.2228 /xm 2000 Jan 11 11:43 3600 

2.2300 Atm 2000 Jan 12 13:26 3600 

2.3129 yum 2000 Jan 11 13:24 4800 



Table 2: Magnetic Field Measurements 









n. 


^eff 


logs 


V sin i 










Bot, (kG) 




Be, (kG) 


Object 


Spectral Type 


(Me) 




(K) 




(km s-') 


(days) 


n 


(10^° ergs s-i; 


1 Ml Fit (x?) 


M2 Fit ixl) 


M3 Fit ixl) 




TWA 5A£ 


Ml.o 






3700 


3.70 


38.0 








4.2 (6.59) 


4.2 (6.49) 


4.9 (6.18) 


1.2 


TWA 7 


M3 


0.92 


1.89 


3300 


3.8-5 


9.0 


5.05 


28 


0.875 


1.6 (1.94) 


2.0 (1.73) 


2.3 (1.88) 


1.3 


TWA 8A 


M2 


0.60 


1.29 


3400 


4.00 


4.0 


4.65 


16 


0.816 


2.3 (1.37) 


2.7 (1.20) 


3.3 (1.25) 


1.4 


TWA 9A 


K7 


0.89 


0.88 


4000 


4..50 


9.0 


5.10 


90 


0.679 


2.4 (2.44) 


2.9 (2.03) 


3.5 (1.74) 


2.1 


TWA 9B 


Ml 


0.30 


0.91 


3400 


4.00 


10.0 


3.98 


60 




2.8 (1.14) 


3.3 (1.06) 


3.1 (1.18) 


1.4 


TW Hya^ 


K7 


2.31 


0.9fi 


4126 


4.84 


.5.8 


2.80 




1.400 


2.2 


2.7 


2.7 


2.1 



"Taken from Lawson fc Crause (2005) . 

''Inclination angles are inferred from rotaion period, stellar 
radius and v sinz. 

'^Observed X-ray luminosities are derived from Webb et al. 
(1999). 

'^Binary or possible triple system. See discussion in Sec. 3.1. 
'^From Yang et al. (2005) . 
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Fig. 1. — Best fit of infrared spectra of TWA 5Aab system using the two-magnetic- 
component model M2. Top two panels: magnetically sensitive Ti I lines; Bottom panel: 
magnetically insensitive CO lines. (Histogram: data; dash-dotted line: fit without magnetic 
field; smooth line: fit with magnetic broadening.) 
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2.220 2.221 2.222 2.223 2.224 
Wavelength()im) 



Fig. 2. — Spectra of TWA 7 along with several spectral type standards in the 2.2228 fim 
setting. The spectra are all vertically offset for better display. Note that bot h AD Leo and G J 



729 are known to have strong magnetic fields ( Johns-Krull fc Valenti 



measu rements for AD Leo by 



Saar fc Linskvl (|l985l ). and GJ 729 by 



2000). See also earlier 



Johns-Krull fc Valenti 



(119961 ). 
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Fig. 3. — Best fit of infrared spectra of TWA 7 using the two-magnetic-component model M2. 
Top two panels: magnetically sensitive Ti I lines; Bottom panel: magnetically insensitive 
CO lines. (Histogram: data; dash-dotted line: fit without magnetic field; smooth line: fit 
with magnetic broadening.) 
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Fig. 4. — Best fit of infrared spectra of TWA 8A using the two-magnetic-component model 
M2. Top two panels: magnetically sensitive Ti I lines; Bottom panel: magnetically insensitive 
CO lines. (Histogram: data; dash-dotted line: fit without magnetic field; smooth line: fit 
with magnetic broadening.) 
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Fig. 5. — Best fit of infrared spectra of TWA 9A using the two- magnetic-component model 
M2. Top two panels: magnetically sensitive Ti I lines; Bottom panel: magnetically insensitive 
CO lines. (Histogram: data; dash-dotted line: fit without magnetic field; smooth line: fit 
with magnetic broadening.) 
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Fig. 6. — Best fit of infrared spectra of TWA 9B using the two-magnetic-component model 
M2. Top panel: magnetically sensitive Ti I lines; Bottom panel: magnetically insensitive 
CO lines. (Histogram: data; dash-dotted line: fit without magnetic field; smooth line: fit 
with magnetic broadening.) 
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Fig. 7. — Observed X-ray luminosity plotted against the predicted X-ray luminosity for the 
TWA stars. The solid line is the line of equality, and the dashed line is where the predicted 
value is 10 times the observed one. The filled circles re present the TWA s tars from this 
paper, with the observed X-ray luminosity derived from 



diamonds represent the Taurus & Auriga stars from 



Webb et al. 



Johns-Krulll (120071). 
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